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Abstract 

We try to understand the quasi-periodic oscillations (QPOs) in low-mass neutron-star and black-hole 
X-ray binaries by a resonance model in warped disks with precession. Our main concern is high-frequency 
QPOs, hectohcrtz QPOs, and horizontal-branch QPOs in the z sources and the atoll sources, and the 
correponding QPOs in black-hole X-ray binaries. Our resonance model can qualitatively, but systematically, 
explain these QPOs by regarding hectohcrtz QPOs as a precession of warp. 
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1. Introduction 

Quasi-periodic oscillations (QPOs) have been observed 
in many low-mass X-ray binaries. They give important 
clues to understand disk structures as well as to evaluate 
the mass and spin of the central neutron stars or black 
holes. Although the mechanism of the QPOs is still under 
debate, recent observations suggest that they can be at- 
tributed to disk oscillations. Observations further suggest 
that some resonant processes are involved in the mech- 
anism of the QPOs, and many disk oscillation models 
in this direction have been proposed since Abramowicz 
and Kluzniak (2001), e.g., Lamb and Miller (2003) and 
Kluzniak et al. (2004). In previous papers (Kato 2003, 
2004a, 2004b) we have proposed that the QPOs are disk 
oscillations excited on a warped disk by a resonant pro- 
cess. The purpose of this paper is to examine how much 
the warp model is compatible with observations. 

We consider a relativistic warped disk. On the disk we 
superpose disk oscillations. Then, some of these disk os- 
cillations have resonant interactions with the disk at par- 
ticular radii through non-linear coupling with the warp 
(Kato 2003, 2004a, 2004b). When the warp has no pre- 
cession (this is the case considered in the above papers) 
the resonances occur at radii where one of the relations of 
K = (V2— l)f2, K = fl/2, or k= (\/3 — is satisfied, where 
k{t) is the epicyclic frequency and f2(r) is the angular fre- 
quency of disk rotation, r being the radius from the disk 
center. In the case of relativistic Kcplerian disks, these 
radii are, in turn, r = 3.62r g , r = 4.0r g , and r = 6.46r g 
(Kato 2003, 2004a), where r g is the Schwarzschild ra- 
dius, defined by r g = 2GM/c 2 , M being the mass of the 
central object. Kato (2004b) subsequently showed that 
among the resonant oscillations mentioned above, those 
at k = f2/2 are excited spontaneously by the resonance 
process, itself. Kato (2004b) further showed that the high- 
frequency QPOs in black-hole X-ray binaries, which are 
usually a pair and have a frequency ratio close to 2 : 3, 



can be explained by this warped model. 

In this paper we extend the warped disk model to the 
case where the warp has precession. Wc demonstrate that 
by this extension the main important characteristics of 
QPOs in X-ray binaries (neutron-star and black-hole X- 
ray binaries) can be qualitatively explained. Among three 
types of resonances, which tend, in the limit of no preces- 
sion, to K= (v2 — K = ^/2, or k = (V3 — l)fl, we focus 
our attention in this paper on the middle one, since the 
resonances in this case spontaneously excite oscillations, 
as mentioned above. 
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Details of the resonance process on warped disks are 
presented by Kato (2003, 2004a, 2004b) in the case where 
the warp has no precession. The essence of the resonance 
processes is the same, even when a warp has precession. 
We thus present here only an outline. An overview of our 
non-linear resonance model in the case where the warp 
has precession is sketched in figure 1 of Kato (2004c). 

We consider geometrically thin disks rotating with an- 
gular velocity Q(r). The epicyclic frequency on the disk 
is denoted by n{r). The oscillations on geometrically thin 
disks are generally classified into g-mode and p-mode os- 
cillations (see, e.g., Kato et al. 1998; Kato 2001). In 
simplified disks the oscillations are further classified by 
the set of (m, n), where m =(0, 1, 2...) is the number of 
nodes in the azimuthal direction, and n =(0, 1, 2...) is a 
number related to nodes in the vertical direction. That is, 
n represents the number of nodes that u r (the radial com- 
ponent of velocity associated with oscillations) has in the 
vertical direction. It is noted, however, that u z (the ver- 
tical component of velocity associated with oscillations) 
has (n — 1) nodes in the vertical direction, and u z = in 
the case of n = 0. 

A warp is a global deformation of disks with m = n = 1. 
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The warp is assumed to have a precession whose angular 
frequency is u> p , On a disk deformed by the warp we su- 
perpose g-modc oscillations with arbitrary m and n. A 
g-modc oscillation with frequency ui and (m, n) has a rel- 
atively large amplitude, global pattern only around the 
radius where 

(w - mfi) 2 - n 2 = (1) 

is satisfied. This can be understood if the dispersion rela- 
tion for local perturbations is considered (e.g., Kato et al. 
1998; Kato 2001). That is, the region of (uj~mti) 2 ~n 2 > 
is a evanescent region of the oscillations. In the region 
where (u — mfl) is smaller than k 2 , on the other hand, 
the oscillations have very short wavelengths in the radial 
direction in geometrically thin disks. 

A non-linear interaction of this g-mode oscillation with 
the warp produces an oscillation with lu±lo p , in and h, 
where m = m ± 1 and h = n ± 1 (these oscillations are 
called hereafter intermediate oscillations). These inter- 
mediate oscillations resonantly interact with the disk at 
the radius where the dispersion relation for these interme- 
diate oscillations is satisfied [see Kato (2004b) for detailed 
discussions]. There are two types of resonances. One is 
resonances that occur through motions in the vertical di- 
rection, and the other is those through motions in the ra- 
dial direction (see Kato 2004a, referred to Paper I). Here, 
we are interested in resonances in the horizontal direction. 
The horizontal resonances occur around the radius where 

(uiwp-mO) 2 -k 2 -0 (2) 

is satisfied. Combining equations (1) and (2), we find that 
the resonances occur at the radius of K= Vt/2±uj p /2 (cf. 
Paper I). After this resonance the intermediate oscillations 
feedback to the original oscillations, amplifying or damp- 
ening the original oscillations (Kato 2004b). Hereafter, we 
consider the case of k = 0/2 + u) p /2. 

A detailed examination shows that when fh of the inter- 
mediate oscillation is m — 1, i.e., m = m — 1, the oscillations 
that resonantly interact with the disk at k = (f2 + w p )/2 
are those satisfying u> — mfl — k at the resonant radius. 
On the other hand, the oscillations that resonantly inter- 
act with the disk at k = (Q + tu p )/2 are those satisfying 
cu = mil + k there, when m = m + 1 (see Paper I) . 

Here, we consider non-axisymmctric oscillations. 
Among them we are particularly interested in oscillations 
of a small number of m. Typical ones are cu = Q— k (m = 1, 
m = 0), lo = Q + k (m = 1, in = 2), and u> = 2£! — k (m = 2, 
m = 1). As shown below, we identify these oscillations, 
in turn, to the horizontal branch QPOs, upper-frequency 
kHz QOPs, and lower-frequency kHz QPOs in the case 
of z sources. Considering this, we introduce the following 
notations: 

luh = Q + k, cjl = 2f2 — k, whbo = O - k. (3) 

3. Precession of Warps 

We consider a relativistic Keplcrian disk with the 
Schwarzschild metric, and examine the radii where the 
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resonance condition, k = (Q + lu p )/2, is satisfied as func- 
tions of uj p . The condition is satisfied at two different radii 
when Up > 0. (In the case of lo p = 0, we have only one ra- 
dius, i.e., 4.0 r g . The other one is oo.) As oj p increases, the 
inner radius becomes larger than 4.0 r g , while the other 
one decreases from infinity. At a certain critical value of 
uj p , both radii coincide and above the critical value of u> p , 
there is no solution of k = (f2 + ui p )/2. The results are 
shown in figure 1. The unit of the abscissa is lo p (M /Mq). 
The critical value of u p is ~ 325 Hz when M/Mq = 1.0, 
while ~ 162 Hz when M/Mq = 2.0. The branch of larger 
value of r on the r-uj p plane is hereafter called the upper 
branch and that of the lower one the lower branch. The 
frequencies (uh, wl, ^hbo) at radii of the lower branch 
decrease as uj p increases, since the resonance radius moves 
outward. On the other hand, the frequencies at radii of 
the upper branch increase as u p increases. To compare 
with observations, the wl-wh, <^hbo -w h and w p -wh rela- 
tions are shown in figure 2, which is free from precession. 
For a comparison, the straight line of the wh-wh rela- 
tion is also added. Values of uj p (M/Mq) are shown, for 
convenience, at some points on the u p -uh curve.. The 
trancated curve on the upper-right corner is the first har- 
monic of whbOj i-C, it represents the 2ajnBO~ UJ n relation. 
The curve extends to the left, but is cutted in order to 
avoid complexity of the figure. It is useful to compare 
this figure with figure 2.9 of van der Klis (2004). We take 
the standpoint that the variation of the precession is the 
cause of time variations of QPO frequencies in a single 
object. 

As the precession frequency changes, the frequencies 
wh, ojl, and whbo vary along the curves in figure 2. 
Observations show that the changes of the horizontal 
branch QPOs and lower-frequency kHz QPOs are corre- 
lated so that the former frequencies are ~ 0.08-times the 
latter ones (Psaltis et al. 1999). Hence, for a compar- 
ison, the curve of the 0.08 wl-wh relation is shown in 
figure 2. Figure 2 shows that the curve of the 0.08 wl- 
ojh relation crosses the curve of the whbo _w h relation 
at ojh ~ 320(M/Mq) _1 Hz, which corresponds to lo p ~ 
115(M/M Q )- 1 Hz. 

It is noted that the curves shown in figure 2 are mass- 
independent, since the axes are normalized by M/Mq. 
That is, the results hold in a wide range of frequency by 
changing the mass. 

4. Summary and Discussion 

In luminous neutron-star low-mass X-ray binaries (the 
z sources), we typically have four distinct types of QPOs. 
These are the ~ 5-20 Hz normal branch oscillation (NBO), 
the 15-60 Hz horizontal branch oscillation (HBO), and 
the - 200-1200 Hz kilohertz QPOs that typically occur 
in pairs. The QPOs in the atoll sources (less luminous 
neutron-star low-mass X-ray binaries) can also be classi- 
fied into similar types of oscillations. In addition, in the 
atoll sources the hectohertz QPOs are observed in the fre- 
quency range of 100 Hz-200 Hz. In the z sources, however, 
the presence of hectohertz QPOs is uncertain, or ambigu- 
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ous. In black-hole X-ray binaries, we have high-frequency 
QPOs in the range of 100 Hz to 450 Hz, which usually 
appear in pairs. A comprehensive review on QPOs is pre- 
sented by van der Klis (2004). 

One of important characteristics of QPOs in the z 
sources is the presence of a strong correlation between 
kilohcrtz QPOs and HBOs. That is, the frequency of the 
lower kHz QPO and that of HBO are correlated in each 
object so that the former is ~ 0.08-times the latter. This 
correlation extends from neutron-star to black-hole X-ray 
binaries in nearly three orders of magnitude in frequency 
(Psaltis et al. 1999, see figure 2 of their paper). In our 
resonance model this correlation can be explained if varia- 
tion of precession frequency occurs around 115(M/M Q ) _1 
Hz, say, 70(M/M Q )~ 1 Hz - 17O(M/M )- 1 Hz. A ques- 
tion is whether such a precession has been observed. In 
relation to this issue, it it interesting to note that the 
observed hectohertz QPOs are roughly in the frequency 
range mentioned above. This suggests that the hecto- 
hertz QPOs might be a manifestation of the warp. One 
of the characteristics of the hectohertz QPOs is that their 
frequencies are nearly constant (e.g., van der Klis 2004). 
In our resonance model, the relation of whbo ~ 0.08 Wl is 
realized over in a wide range of wh without much chang- 
ing the value of u> p , as shown in figure 2. This is con- 
sistent with the idea that the warp represents the hec- 
tohertz QPOs. The precession of disks in X-ray binaries 
is theoretically expected, since the radiation force from 
central star gives torques on warped disks (Pringle 1996; 
Maloney et al. 1996). It is not clear, however, whether 
such a high-frequency precession as required here is gen- 
erally expected. We suppose that the frequency of the 
precession is related to rotation of the central star. A 
non-axisymmetric pattern on a rotating stellar surface will 
give rise to a precession of the disk through the effects of 
a radiative force or a magnetic field. 

Another important characteristic of QPOs in neutron- 
star X-ray binaries is that the frequency ratio of the pair 
kHz QPOs is close to 3 : 2, but changes with time so that 
the ratio decreases with an increase of the frequency. This 
observational trend is realized in our model. In our model 
the ratio wh/wl is 1.5 at wh ~ 700(M/M Q ) _1 Hz, and 
becomes larger than 1.5 for smaller cjjj- For larger uh the 
ratio tends to unity as the resonant radii approaches to 
4.0 r g (i.e., lu p = 0). 

Next, let us consider black-hole X-ray binaries. The 
pair of high-frequency QPOs in these objects changes little 
their frequencies, keeping the ratio close to 3 : 2. This is 
different from the case of neutron-star X-ray binaries. Our 
resonance model qualitatively explains this. Considering 
their mass and the observed frequencies of QPOs, we think 
that in the case of black-hole X-ray binaries the observed 
QPOs are those resulting from the lower branch of the r— 
ujp relation (see figure 1). That is, the resonance radius is 
close to 4.0 r g . This case corresponds to the upper- right 
corner of figure 2. In the limit of o; p = 0, the resonance 
occurs at r = 4.0 r g and the frequency ratio of cjl and 
the first harmonic of whbo is just 3 : 2. (It is noted that 
the first harmonic of c^hbo has been also observed in the 



z sources.) These QPOs change little their frequencies 
for a change of precession frequency, since as shown in 
figure 1 the resonance radius remains close to 4.0 r g for a 
large change of u> p . Another explanation of little change of 
frequencies of high-frequency QPOs in black-hole binaries 
is that the precession is really small in the case of black 
holes, since the radiation force from the central object is 
absent. It is noted that in our model the modes of the pair 
QPOs in the black-hole binaries are different from those 
in neutron-star binaries. That is, the pairs in the former 
are wh(=wl) and 2whbOj while those in the latter are <x>h 
and wl- 

We have restricted our attention only to the resonances 
at k = (fi + cJ p )/2, since the analyses of growth rate 
of resonant oscillations suggest that the resonances at 
K = (V2 — l)f2 + w p and k = (v3 — l)f2 + cu p are not spon- 
taneously excited (Kato 2004b). However, examinations 
of these cases are worthwhile. As an example, some re- 
sults in the case of k = (\/2-l)fi + uj p arc shown in fig- 
ure 3. For simplicity, in this figure, only the resonant 
oscillations resulting from the upper branch of the r-w p 
curve are shown. The frequency of precession required 
to obtain cjhbo ~ 0.08 ui^ is around 50(M/M Q ) _1 Hz- 
100(M/M Q )- 1 Hz. 

Our model predicts that some QPOs that are not yet 
observed are present in neutron-star and black-hole X-ray 
binaries. In our model the observed QPOs in neutron-star 
X-ray binaries are the resonance oscillations at resonant 
radii belonging to the upper branch of the r-uj p diagram 
(figure 1). Resonances coming from a resonance radius 
of the lower branch are also expected. They have higher 
frequencies compared with the observed ones. The ob- 
served QPOs in black-hole X-ray binaries, on the other 
hand, are interpreted to be resonant oscillations belong- 
ing to the lower branch of the i — lo p relation. Resonant 
oscillations resulting from radii of the upper branch are 
also expected. These oscillations have lower frequencies 
compared with the observed ones. 

A note added on April 30: 

In the present paper we have considered horizontal res- 
onances of g-mode oscillations. In this model, the pre- 
cession required is retrograde. Furthermore, time varia- 
tion of the precession is required in order to explain time 
variation of QPO frequencies. In the case of vertical reso- 
nances of g-mode oscillations, however, the observed QPO 
frequencies and their time variations can be explained as 
a result of time change of disk structure in the vertical di- 
rection, without appealing to precession. This is discussed 
in a subsequent paper. 
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Fig. 1. The i — uj p relation giving the solution of the reso- 
nance condition k= (Q + lu p )/2. The disk is Kcplerian in the 
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Fig. 2. The uj^-ljh, ojhbo^^H; an d <^p-<^H relations ob- 
tained from the resonance condition k = (Q + w p )/2 for some 
frequency range of ojjj • For a comparison, the line of ojh - 
is shown and the value of u> p along the ui p -ujh curve is shown 
at some points. For a comparison, the 0.08 cjl-^h relation 
is shown. The trancated curve on the upper-right corner rep- 
resents the first hamonic of i^hbOj i- e -> it i s the 2cjhbo — U E 
relation. In order to avoid complexity, the curve is trancated. 
It is noted that the ratio ujj/^l decreases as ljjj increases. 
The ratio is ~ 1.5 for ton in the middle of the figure, and 
becomes 1.0 at the right end of the figure, where ojh is the 
maximum and lu p = 0. At this right end, the ratio i^h(= ^l) 
: 2u;jjbo : ^HBO = 3 : 2 : 1. We regard the oscillations of 
u} p as hectohertz QPOs. 
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Fig. 3. Same as figure 2, except for k = {y2— l)Q + u; p . In 
this figure, for simplicity, only the resonant oscillations that 
occur on the upper branch of the i — u) p plane are shown. 
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